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Various fundamental phenomena of strongly-correlated quantum systems such as high- Tc super- 
conductivity, the fractional quantum-Hall effect, and quark confinement are still awaiting a univer- 
sally accepted explanation. The main obstacle is the computational complexity of solving even the 
most simplified theoretical models that are designed to capture the relevant quantum correlations 
of the many-body system of interest. In his seminal 1982 paper [I], Richard Feynman suggested 
that such models might be solved by "simulation" with a new type of computer whose constituent 
parts are effectively governed by a desired quantum many-body dynamics. Measurements on this 
engineered machine, now known as a "quantum simulator," would reveal some unknown or difficult 
to compute properties of a model of interest. We argue that a useful quantum simulator must 
satisfy four conditions: relevance, controllability, reliability, and efficiency. We review the current 
state of the art of digital and analog quantum simulators. Whereas so far the majority of the focus, 
both theoretically and experimentally, has been on controllability of relevant models, we emphasize 
here the need for a careful analysis of reliability and efficiency in the presence of imperfections. We 
discuss how disorder and noise can impact these conditions, and illustrate our concerns with novel 
numerical simulations of a paradigmatic example: a disordered quantum spin chain governed by 
the Ising model in a transverse magnetic field. We find that disorder can decrease the reliability 
of an analog quantum simulator of this model, although large errors in local observables are intro- 
duced only for strong levels of disorder. We conclude that the answer to the question "Can we trust 
quantum simulators?" is... to some extent. 



I. INTRODUCTION 

In his 1982 foundational article [T], Richard Feynman 
suggested that the complexities of quantum many-body 
physics might be computed by "simulation." By design- 
ing a well-controlled system from the bottom up, one 
could create a computer whose constituent parts are 
governed by quantum dynamics generated by a desired 
Hamiltonian. Measuring the properties of this nano- 
engineered system thus reveals some unknown or difficult 
to compute properties of a quantum many-body model, 
such as the nature of quantum phase diagrams. Feyn- 
man's machine is now known as a "quantum simulator" 
(QS). 

Fueled by the prospect of solving a broad range of 
long-standing problems in strongly-correlated systems, 
the tools to design, build, and implement QSs [IH3] have 
rapidly developed and are now reaching very sophisti- 
cated levels [4]. Researchers are making breakthrough 
advances in quantum control of a variety of systems, 
including ultracold atoms and molecules ([5til7j. for re- 
views see also [ISfEO] ). ions ( |21ff27] . for recent reviews 
see |20 1 l28H30j ) . photons ( |31ff36) . for a recent review see 
[37]), circuit quantum electrodynamics (CQED) and po- 
laritons ( |38fH^ . for a recent review see ^43]), artificial 
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lattices in solid state |44| . nuclear magnetic resonance 
(NMR) systems [i51 - l5U] . and superconducting qubits (for 
reviews see |51ff54] . for the current state of art [55| [56] 
and references therein). For a general overview see also 
|57| . At the current pace, it is expected that we will 
soon reach the ability to finely control many-body sys- 
tems whose description is outside the reach of a classical 
computer. For example, modeling interesting physics as- 
sociated with a quantum system involving 50 spin-1/2 
particles - whose general description requires 2^" w 10^^ 
complex numbers - is out of the reach of current classical 
supercomputers, but perhaps within the grasp of a QS. 

In a field brimming with excitement, it is important to 
critically examine such high expectations. Real-world im- 
plementations of a quantum simulation will always face 
experimental imperfections, such as noise due to finite 
precision instruments and interactions with the environ- 
ment. Feynman's QS is often considered as a fundamen- 
tally analog device, in the sense that all operations are 
carried out continuously. However, errors in an analog 
device (also continuous, like temperature in the initial 
state, or the signal-to-noise ratio of measurement) can 
propagate and multiply uncontrollably [581. Indeed, Lan- 
dauer, a father of the studies of the physics of informa- 
tion, questioned whether quantum coherence was truly 
a powerful resource for computation because it required 
a continuum of possible superposition states that were 
"analog" in nature |59| . 

This contrasts with the operation of a universal digi- 
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tal quantum computer as envisioned by David Deutsch. 
in which all operations are digitized into a finite set of 
logic gates and measurements [BD] [BT]. The invention 
of quantum-error-correcting codes showed that a quan- 
tum computer is in some sense both analog and digital. 
Through a discrete set of unitary transformations, we can 
get arbitrarily close to any superposition, and imperfec- 
tions can always be projected on a discrete set and thus 
can be corrected [52] . When such a digital quantum sim- 
ulation operates with fault-tolerant quantum error cor- 
rection [63], we can trust its output to a known finite 
precision. 

Universal digital quantum computers may serve as dig- 
ital QSs (DQSs) that mimic dynamics of some quantum 
many-body system of interest. Despite the fact that in 
such a case error correction and fault tolerance is guaran- 
teed, the question of efficiency of such a device is highly 
non-trivial. The number of resources needed for precise 
simulation of continuous-time dynamics of a many-body 
system by stroboscopic digital applications of local gates 
might be enormous One can also consider DQSs 

that are experimental systems that have at their disposal 
only a limited, non-universal set of gates. In such a sit- 
uation, the error correction and fault tolerance are not 
guaranteed and the question of efficiency is even more 
pertinent. 

This raises the central problem of this key issue arti- 
cle: can we trust the results obtained with a real- world 
analog or digital QS, and under what conditions are they 
reliable to a known degree of uncertainty? Although our 
main discussion concentrates on analog QSs (AQSs) , the 
article reports also on the state of art of DQSs. It is orga- 
nized as follows. Section II develops the general concept 
of QSs in the spirit of the DiVincenzo criteria for quan- 
tum computing [63]. Here we present one of the main 
results of this article: a definition of the QS based on 
four properties that a QS should have: relevance, con- 
trollability, reliability, and efficiency. 

Section III is devoted exclusively to DQSs. It contains 
several subsections in which we review various proposals 
for DQSs, classify them, and discuss the present state 
of knowledge concerning their controllability, reliability, 
and efficiency. Section IV is organized similarly, but fo- 
cused on AQSs. 

Section V is perhaps the most important one from the 
conceptual point of view. Here, we formulate specific 
proposals to investigate the robustness of AQS and how 
to extend standard methods of validation and certifica- 
tion of AQS. We illustrate these considerations in Sec- 
tions VI- VIII with calculations for a paradigmatic model 
that are not published elsewhere. Section VI describes 
the investigated model, and the results concerning stat- 
ics and dynamics are presented in Sections VII and VIII, 
respectively. We conclude in Section IX. The paper con- 
tains also an Appendix describing technical details of the 
methods used. 

Very recently. Nature Physics has published a focus 
issue with five articles devoted to QSs: a short article by 



J.I. Cirac and P. Zoller [BB], introducing the subject, and 
four longer reviews on ultracold atoms |19j, ions |30j, 
CQED [67J, and photons [37]. Our key issue article is 
complementary to these, in the sense that it addresses 
general and universal problems of validation of quantum 
simulations, their robustness, reliability, and efficiency - 
problems that pertain to all kinds of QSs. 



II. QUANTUM SIMULATORS 

Before proceeding, we must establish a clear definition 
of a QS. We consider here a QS to be a device which, 
when measured, reveals features of an ideal mathemati- 
cal model, e.g., the phase diagram for the Bose-Hubbard 
model on a specified lattice with specified interactions. 
This contrasts with, and is less demanding than, a full 
simulation of a real material, since typically a mathemat- 
ical model attempts to capture only the most relevant 
properties of the real material. For example, the super- 
conducting properties of a cuprate might be shared, in 
part, by a Fermi-Hubbard model |68. i69j . A QS may be 
a special purpose device that simulates a limited class 
of models, e.g, the Bose-Hubbard model simulated by 
atom transport in an optical lattice [70, 71^, or a univer- 
sal machine that is capable, in principle, of simulating 
any Hamiltonian on a finite-dimensional Hilbert space. 

Based on this, we formulate the following "working" 
definition of a quantum simulator in the spirit of the Di- 
Vincenzo criteria for quantum computing |65| , providing 
some more detailed explanations below (see also [3J): 

Definition 

A QS is an experimental system that mimics a simple 
model, or a family of simple models of condensed mat- 
ter (or high-energy physics, or quantum chemistry, . . . ). 
A quantum simulator should fulfill the following four re- 
quirements: 

• (a) Relevance: The simulated models should be 
of some relevance for applications and/or our un- 
derstanding of challenges in the areas of physics 
mentioned above. 

• (b) Controllability: A QS should allow for broad 
control of the parameters of the simulated model, 
and for control of preparation, initialization, ma- 
nipulation, evolution, and detection of the relevant 
observables of the system. 

• (c) Reliability: Within some prescribed error, 
one should be assured that the observed physics of 
the QS corresponds faithfully to that of the ideal 
model whose properties we seek to understand. 

• (d) Efficiency: The QS should solve problems 
more efficiently than is practically possible on a 
classical computer. 
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digital analog 

universal non-univ. open universal non-universal open 



Realizations 



trapped ions, 
ultracold neutral 
or Rydberg atoms, 
circuit QED, super- 
conducting qubits, . 



same as 

universal 

digital 



as universal 
digital (especially 
trapped ions, 
ultracold neutral 
or Rydberg atoms) 



many (trapped ions, 
ultracold atoms, pho- 
tonic and polariton 
systems, artificial 
solid-state lattices, . . 



same as 
non- 
universal 
analog 



Control full (long-range inter- partial partial 

actions difficult ?) 



full partial (but long-range partial 

interactions "easy") 



Error 

correction 

(EC) 



with exponential 
overhead 

(Trotterization issues) 



not 

guaranteed 



not 

guaranteed 



no 

standard 
EC 



no 

standard 
EC 



no 

standard 
EC 



Reliability full 



not 

guaranteed 



not 

guaranteed 



? (partial 
validation 
schemes 
available) 



? (partial 
validation 
schemes 
available) 



Efficiency 



efficient without EC 
(for general class of 
models); much less 
efficient with EC 
(Trotterization issues) 



at least as 
universal 
digital, but 
may not be 
provable 



can be better 
than universal 
digital 



Table I. In this table, we characterize the different classes of quantum simulators (digital vs. analog, universal vs. non-universal, 
Hamiltonian vs. open) focussing especially on the requirements (b) to (d). (Since the relevance (a) depends on the concrete 
model simulated, we do not list it here.) Note specifically that little is known about reliability and efficiency of AQSs (question 
marks). Detailed descriptions are provided in Sections III and IV. 



In Table |T] we summarize to which extent existing 
experimental proposals fulfill these requirements. We 
will characterize the different types of QSs in more detail 
in the following two sections, but before that, we would 
like to make some general comments. 

Comments ad a) We should demand that the mim- 
icked models are not purely of academic interest but that 
they rather describe some interesting physical systems 
and solve open problems. This means also that the 
simulated models should be computationally very hard 
for classical computers (see also requirement (d)). 

Comments ad b) and c) Regarding control over mea- 
surable observables, one should stress that very often the 
amount of output information required from quantum 
simulators might be significantly smaller than one could 
demand from a universal quantum computer. Quantum 
simulators should provide information about phase dia- 
grams, correlation functions, order parameters, perhaps 
even critical exponents or nonlocal hidden order param- 
eters. But a common assumption is that these quanti- 
ties are more robust than what is required for a univer- 
sal quantum computer, which typically relies on much 



higher-order correlation functions than a QS. 

Regarding control over model parameters, it is in 
particular desirable to be able to set the parameters in 
a regime where the model becomes tractable by clas- 
sical simulations, because this provides an elementary 
instance of validating the QS. Furthermore, one of the 
main results of this paper is the proposal and analysis of 
an even more sophisticated manner of validation, namely 
the checking of the sensitivity of the quantum simulation 
with respect to addition of noise and /or disorder. Such 
a validation is only possible with sufficient control over 
the system. Note, however, that there are other possi- 
bilities of checking the results, as pointed out to us by 
Z. Hadzibabic |72]. Namely, sometimes it is impossible 
to simulate the system classically, but it might still be 
possible by classical means to test the sensibleness of the 
quantum-simulation results. For instance, the measured 
ground-state energy should fulfill all known bounds, 
such as variational ones, and others. 

Comments ad d) The notion "computationally very 
hard for classical computers" may have several meanings: 
i) an efficient (scalable, with polynomial growth in re- 
sources as a function of problem size) classical algorithm 
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to simulate the model might not exist, or might not be 
known; ii) the efficient scalable algorithm is known, but 
the required size of the simulated model is too large to 
be simulated under reasonable time and memory restric- 
tions. The latter situation, in fact, starts to occur with 
the classical simulations of the Bose- or Fermi-Hubbard 
models [TIT , in contrast to their experimental quantum 
simulators. However, there might be exceptions to the 
general rules. For instance, it is desirable to realize QSs 
to simulate and to observe novel phenomena that so far 
are only theoretically predicted, even though it might 
be possible to simulate these phenomena efficiently with 
present computers. Simulating and actually observing 
in the lab is more than just simulating abstractly on a 
classical computer. 

Comments ad c) and d) The requirements of relia- 
bility and efficiency are interrelated. In fact, we could 
try to improve the precision of a QS by averaging more 
experiments, but in hypersensitive regimes (like close 
to quantum phase transitions) the necessary number 
of repetitions can grow rapidly, bringing the overall 
efficiency of the QS down to the level of classical com- 
puters. A connection between (c) and (d) could also be 
relevant for the popular cross-validation approach |74| . 
There, one compares the results of two different physical 
realizations performing a QS of the same model, and 
hopes to find universal features which then would be 
ascribed to the simulated model. It may be, however, 
that the universal features shared by multiple platforms 
are robust only because they could have been predicted 
efficiently with some classical algorithm. 

With our working definition at hand, one could ask - 
What should a QS simulate? An important set of tasks 
include: 

1. Statics of the mimicked system at zero tempera- 
ture; this implies ground-state simulation and its 
properties. 

2. Statics at thermal equilibrium, i.e., Hamiltonian 
dynamics at low energies or thermodynamics at 
non-zero, typically low, temperatures. 

3. Continuous-time dynamics of the system, in partic- 
ular Hamiltonian dynamics out of equilibrium. 

4. Dissipative or open-system continuous-time dy- 
namics. 

To understand which of these are most relevant, we 
now discuss shortly which systems can be simulated effi- 
ciently classically and which systems are classically com- 
putationally hard. Classical simulations of quantum sys- 
tems are currently performed using one of the following 
numerical methods [3J: 

• Quantum Monte Carlo (QMC) 

• Systematic perturbation theory 



• Exact diagonalizations 

• Variational methods (mean field methods, density- 
functional theory (DFT), dynamical mean-field 
theory (DMFT), tensor-network states (TNS), 
density-matrix renormalization group (DMRG), 
tree tensor network states (TTN), multiscale 
entanglement-renormalization ansatz (MERA) , 
projected entangled-pairs states (PEPS), ...) 

Each of these methods has its limitations. Let us first 
focus on points 1) and 2) of the previous list of possible 
QS tasks. In these cases, QMC works for various large 
systems, but fails for Fermi or frustrated systems due 
to the famous sign problem |7S]. Perturbation theory 
works only if there exists a small expansion parameter 
|76| . Exact diagonalization works only for rather small 
systems (T^. In the case of ID systems, DMRG, MERA 
and TTN techniques scale favorably and can, in principle, 
treat very large systems [TTfTf^] . In 2D the situation is 
more complex - similar to exact diagonalization, DMRG 
and TTN work only for reasonably small systems |80ff5^ , 
whereas 2D tensor-network methods (PEPS, MERA) in 
principle work for arbitrarily big systems (bosonic, and 
even fermionic [83 or frustrated [84 J but are biased to- 
wards slightly entangled states. Mean field [85 , DFT 
[HSl IH7], or DMFT [88 , finally, have other limitations, 
e.g., they are essentially designed for weakly-correlated 
systems. 

Which are then the models that are computationally 
hard for points 1) and 2) in the previous task list? Gen- 
erally speaking, computationally hard are those "strongly 
entangled" models in more than ID such as 

• Fermionic models, with paradigmatic examples be- 
ing the Fermi-Hubbard or t — J models for spin 1/2 
fermions [68J. 

• Frustrated models, with paradigmatic examples be- 
ing antiferromagnetic Heisenberg or XY models on 
a kagome or anisotropic triangular lattice |89) . 

• Disordered models, with paradigmatic models be- 
ing quantum, or even classical spin glasses |90] . 

When we move to points 3) and 4) of the task list, i.e., 
studying dynamics, one can safely state that 

• Quantum dynamics on a long time scale is generi- 
cally computationally hard. 

The latter statement implies that while it might be pos- 
sible to simulate with classical computers short-time dy- 
namics in a restricted class of ID models, such attempts 
will nearly always fail at longer time scales. Indeed, this 
fact is related to correlation and entanglement spreading 
according to the Lieb-Robinson theorem that states that, 
after a sufficiently large time, states can become strongly 
entangled ([HHMI, see also [57]). 

In the following two sections, we will explore in more 
detail the state of the art concerning the four require- 
ments (a-d) of our definition, first for digital, then for 
analog QSs. 
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III. DIGITAL QUANTUM SIMULATORS (DQS) 

In this section, we classify DQSs, discuss their general 
properties, various protocols for implementing such de- 
vices, and summarize state-of-art knowledge concerning 
their controllability, reliability, and efficiency. 

A. Universal Digital Quantum Simulators (UDQS) 

While the concept of QSs should be traced back to 
prophecies of Feynman [T] , the ideas were made concrete 
by Lloyd who showed that any "local" many-body uni- 
tary evolution governed by a "local" Hamiltonian could 
be implemented by the control afforded by a universal 
digital quantum computer |98j . For this reason, in the 
following we will term Lloyd's DQS a "universal DQS" 
(UDQS). 

Lloyd's UDQS is in fact a universal quantum computer, 
whose task is to simulate the unitary time-evolution op- 
erator of a certain quantum system described by a phys- 
ical Hamiltonian, which can then be employed to extract 
quanties like energy gaps and ground-state properties. 
This is done by appropriate subsequent stroboscopic ap- 
plications of various quantum gates that mimic the action 
of a global unitary continuous time evolution operator of 
the system. The mathematical basis for such a digit al- 
ization is given by the Trotter-Suzuki formula. In order 
to realize Lloyd's UDQS in a laboratory, the experimen- 
talist has to have to his/her disposal a universal set of 
unitary quantum gates [HH!- Let us list below some pos- 
sible realizations and properties of UDQSs: 

• Realizations: While implementation of a fully- 
functioning large-scale digital quantum computer is 
still in development, there are several physical sys- 
tems for which the universal sets of quantum gates 
are available, and for which realization of proof- 
of-principle UDQSs is possible. These systems in- 
clude ultracold ions [30', ultracold trapped atoms 
interacting via cold collisions [19 , or the Rydberg- 
blockade mechanism jlOOl ITOT]. circuit QED [67] . 
superconducting qubits (for reviews see |CTI5i] , for 
the current state of art and references therein). 
For a general overview see also [ST. The first 
concrete proposals for realization of UDQSs where 
given in |102U103| , and the first experiments, per- 
haps, were performed in NMR systems [JSHIZI EDI] ■ 
Using a digital architecture and stroboscopic se- 
quence of gates, the quantum simulation of Ising, 
XY, and XYZ spin models in a transverse field were 
recently demonstrated in a proof-of-principle exper- 
iment with up to six ions |24| . 

• Controllability: In accordance with [5S], a UDQS 
is perfectly controllable, i.e., with the help of a uni- 
versal set of gates sufficient control of the parame- 
ters can be achieved. This control allows for simula- 
tion of practically any local Hamiltonian evolution. 



as well as for preparation, manipulation, and detec- 
tion of relevant states and observables of the system 
in question. Further, Preskill's group has proven 
recently that the scattering amplitudes in the sim- 
ple relativistic quantum field theories can be effi- 
ciently (in polynomial time) simulated by UDQSs 
|105II106| . Note, however, that not much is known 
about the possibility of quantum simulation of sys- 
tems with long-range interactions like Coulomb or 
dipole-dipole interactions using UDQSs. 

• Error correction: A UDQS is the only DQS 
which has guaranteed access to error correction 
and fault tolerance |1071 1108] (for the first proof- 
of-principle experiments see |56l I109flll3| l. 

• Efficiency: So far, the community has mostly fo- 
cused on developing requirement (b) for suitable 
relevant models, both theoretically and experimen- 
tally. The conditions (c) and (d) have received 
considerably less attention, especially their inter- 
relation. Most work is focused on efficiency in the 
absence of errors. Lloyd showed that a Trotter- 
Suzuki decomposition of a time-evolution operator 
is efficient in that each logic gate acts on a scal- 
able Hilbert space associated with a small subset 
of qubits and the total number of gates N scales 
polynomially, N ~ t^/e, where t is the time of evo- 
lution to be simulated and e is the error in the re- 
sult |98]. Aharanov and Ta-Shma showed that a 
UDQS is efficient when the Hamiltonian is "sparse," 
i.e., the number of nonzero entries in any row is at 
most poly(log(D)), where D is the dimension of 
the many-body Hilbert space [114' . In the absence 
of errors, the computational complexity of such a 
simulation has been well studied [1151 1116] . 

• Reliability: In the presence of errors, however, en- 
suring reliability to a desired precision has profound 
implications for efficiency even in a digital simula- 
tor on a fault-tolerant quantum computer [M], I117j . 
In the digital approach with a finite universal gate 
set, one applies error-correction schemes that can 
make the whole computation fault-tolerant when 
the error per operation is below a certain thresh- 
old - thus digital simulators fulfill the reliability 
requirement (c). The Trotter expansion, however, 
can scale poorly when error correction is included, 
as emphasized by Brown et al. |64] in studies of 
an implementation of a quantum algorithm to cal- 
culate the low-lying energy gap in pairing Hamil- 
tonians |118| . Because the number of gates in 
the expansion scales as 1/e, in order to achieve 
M = — log2(e) bits of precision, we must Trotter- 
ize the unitary evolution to be simulated into 2^ 
slices. In the presence of errors, each of the time 
slices must be implemented with only a finite set 
of universal gates, according to the Solovay-Kitaev 
theorem [119]; only then can they be implemented 
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fault tolerantly. The result is that a fault-tolerant 
implementation of the Trotter expansion requires 
a number of gates and time to perform the simu- 
lation that grows exponentially with the degree of 
precision required, for a fixed number of particles 
being simulated. Moreover, Brown et al. showed 
that for a small number of qubits where one might 
avoid error correction, analog control errors on the 
logic gates can lead to faulty results, negating re- 
quirement (c), and robust control pulses become 
essential. 

In a similar vein, Clark et al. |120| performed a 
careful analysis of the resources necessary to im- 
plement the Abrams-Lloyd algorithm (121) to cal- 
culate the ground-state energy of the one dimen- 
sional transverse Ising model (TIM) using a state- 
of-the-art fault-tolerant architecture for an ion-trap 
quantum computer. Again, the overhead in the 
number of time steps to fault-tolerantly implement 
the quantum-phase estimation algorithm grows ex- 
ponentially with the degree of precision required. 
They found that for 100 spins, in order to achieve 
6 > 10 bits of precision, at least two levels of con- 
catenated error correction are necessary, requiring 
at least 100 days of run time on the ion-trap quan- 
tum computer; for b > 18, three levels are nec- 
essary, requiring at least 7.5 x 10'^ years! These 
results assume a gate time of 10/iS. To recover the 
100 days limit with only one level of concatenated 
error-correction coding, a gate time of 300ns seems 
necessary (as well as decreasing other parameters 
such as failure probabilities). On the other hand, 
for a fixed precision, the number of resources re- 
quired grows weakly with system size. So, if the er- 
ror probability per gate can be reduced well below 
threshold to achieve the desired precision without 
many layers of concatenated error-correction en- 
coding, then digital quantum simulation will scale 
favorably with the number of particles. 

Let us finally remark that, to increase their effi- 
ciency, digital-quantum-simulation algorithms often com- 
press the number of degrees of freedom that are necessary 
to describe the many-body system, rather than directly 
map the Hilbert space of the system to the Hilbert space 
of the simulator |117| . an approach that has been bor- 
rowed from classical algorithms like MPS or PEPS. Cur- 
rently, there is a new theoretical development towards 
a "hybrid" device, where the ground state of many-body 
Hamiltonians is represented as a PEPS, but implemented 
on a quantum computer. This is efficiently possible when 
the gap between the ground and first excited state scales 
as the inverse of a polynomial in the number of particles. 
Then, one can use the quantum computer to contract 
tensor networks and use that to calculate the expecta- 
tion value of any local variable, such as correlation func- 
tions |122| . In a similar spirit, Temme et al. developed 
a quantum-algorithmic version of the Metropolis Monte- 



Carlo algorithm that allows one to efficiently sample from 
a Gibbs thermal state |123| . Such approaches point to ef- 
ficient DQSs for well defined classes of problems. 



B. Non-Universal Digital Quantum Simulators 
(nUDQS) 

A non-Universal Digital Quantum Simulator (nUDQS) 
is in many aspects similar to a UDQS, except that is it a 
special-purpose quantum computer. Its task is, however, 
the same as that of a UDQS: to simulate continuous-time 
quantum many-body dynamics of a certain quantum sys- 
tem described by a certain physical Hamiltonian. The 
experimentalist who realizes a nUDQS has to his/her dis- 
posal a non-universal set of unitary quantum gates. Let 
us list below some properties and possible realizations of 
such a nUDQS: 

• Realizations: In all systems in which the univer- 
sal sets of quantum gates are available, one can 
also restrict the set of gates and realize a nUDSQ. 
For example, in some of the recent experiments of 
Blatt j25 , only a necessary subset of the available 
set of universal gates was used. All of the systems 
discussed above (atomic, superconducting, etc.) are 
potentially platforms for implementing nUDQSs. A 
seminal example of this approach goes back to the 
so-called "Average Hamiltonian Theory" in nuclear 
magnetic resonance (NMR) |1241 1125] . 

• Controllability: nUDQSs are typically not per- 
fectly controllable, but in most experimental real- 
izations should allow for a wide control of parame- 
ters, which in turn should allow for simulations of 
evolution for wide families of Hamiltonians of in- 
terest. 

• Error correction: For nUDQSs, it is not guaran- 
teed that error correction and fault-tolerant com- 
puting is possible. 

• Efficiency and reliability: All of the above dis- 
cussion concerning UDQSs applies also to nUDQSs. 
But, there are many novel, open problems associ- 
ated specifically with nUDQSs, since, e.g., some- 
times giving up on universality can result in sub- 
stantial efficiency gains. For example, universality 
could be sacrificed in favor of a highly precise and 
fast gate (a simple example is an external ho- 
mogeneous field, which in a UDQS might have to 
be applied as a sequence of one qubit gates). In 
particular, it is possible that for some classes of 
nUDQSs the problems of Trotterization are not as 
severe as in the case of UDQSs [101) . 
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C. Open-System Digital Quantum Simulators 
(OSDQS) 

An Open-System Digital Quantum Simulator (OS- 
DQS) is a completely new concept, in principle very dif- 
ferent from DQSs aimed at Hamiltonian evolutions. OS- 
DQSs are designed to simulate open-system, dissipative 
dynamics described in the simplest situation by a Marko- 
vian Lindblad master equation for the density matrix of a 
many-body system of interest. OSDQSs can be aimed at 
a continuous-time simulation of interesting open-system 
dynamics, or at a designed dissipative dynamics toward 
a stationary state of interest, in particular a pure, highly- 
entangled state [nsHnH]- 

The experimentalist who realizes an OSDQS, in con- 
trast to a UDQS or a nlJDQS, needs to have at his/her 
disposal some non-unitary, dissipative quantum gates, 
which mathematically correspond to Lindblad super- 
operators acting on the density matrix in the master 
equation. This fact opens a plethora of new questions, 
e.g., what are the universal sets of gates for this type of 
evolution. Note that in the case of unitary computing, 
the universal set of gates allows for realization of arbi- 
trary unitary transformations acting on the (pure) state 
of the system. In the case of open-system dynamics, a 
universal set of gates should allow for the realization of an 
arbitrary completely positive map (CPM) acting on the 
density matrix of a system. Moreover, for experimental 
realizations, we require the gates to be local. 

While the conditions for controllability of an open 
quantum systems are under exploration |129] . the ques- 
tion of a universal set of gates in this context remains 
open. A non-trivial reduction (cf. [130, 131 ) of this ques- 
tion to the CPMs that correspond to Markovian evolu- 
tion, is also open. The problem of error correction in 
this context is unsolved as well. All of these comments 
imply that in the area of OSDQSs there are more open 
questions than answers. 

• Realizations: In systems in which the universal 
set of quantum gates is available, one way to real- 
ize dissipative gates is by tracing out ancillas, thus 
allowing to realize an OSDQS. Good testbeds for 
exploring OSDQSs are provided by Rydberg atoms, 
atomic ensembles, NMR {132] , or trapped ions. In 
fact, the first concrete proposals for open system 
DQS concerned Rydberg gates [1001 [M]- The first 
experimental realizations of these ideas, however, 
have been achieved with trapped ions [25;. 

• Controllability: OSDQSs are typically not uni- 
versal since they are not usually controllable in 
the sense of realizing an arbitrary quantum map 
|133| . Nevertheless, many experimental realiza- 
tions should allow for a wide control of parame- 
ters, which in turn should allow for simulations of 
open-system (Markovian) evolutions for wide fam- 
ilies of open systems of interest. As pointed out 
in |127l 1128] , due to the purely dissipative nature 



of the process, this way of doing quantum informa- 
tion processing exhibits some inherent robustness 
and defies some of the DiVincenzo criteria for quan- 
tum computation. In particular, there is a natural 
class of problems that can be solved by open-system 
DQSs or AQSs: the preparation of ground states of 
frustration-free quantum Hamiltonians. 

• Error correction: For OSDQSs, it is not guaran- 
teed that error correction and fault-tolerant com- 
puting is possible in the sense defined above |134| 

• Efficiency and reliability: All of the above dis- 
cussion concerning UDQSs and nlJDQSs applies 
also to OSDQSs. But, due to purely dissipative 
nature of the process, this type of simulation has a 
certain intrinsic robustness and built-in "error cor- 
rection." A clear example is seen in the OSDQS 
implementation of Kitaev's toric code |25l llOlj . 
However, as discussed in Ref. |r01j , errors in the 
gates result in effective heating. Also, the prob- 
lems of Trotterization are not as severe as in the 
case of quantum simulators of Hamiltonian evolu- 
tion. Still, most of these general aspect concerning 
OSDQSs have not yet been investigated systemat- 
ically. The efficiency of OSDQSs for the case of 
frustration-free Hamiltonians depends on the size 
of the gap between the ground state and the ex- 
cited states, or more precisely on the real part of 
the first non-zero eigenvalue of the Lindblad equa- 
tion, which determines the rate of approaching the 
stationary (ground) state. 

Currently, considerable attention has been devoted to 
the problem of existence and uniqueness of the open- 
system preparation of ground states of frustration-free 
Hamiltonians, and in particular, entangled states of in- 
terest. These states are annihilated simultaneously by 
all of the local frustration-free Lindblad superoperators 
entering the master equation. There is little known in 
general about the many-body dissipative dynamics with 
a frustrated set of Lindblad superoperators competing 
with Hamiltonian dynamics. For the first attempts to 
understand these kind of problems in the context of quan- 
tum diffusion-exclusion processes competing with Hamil- 
tonian evolution see, e.g., |135| . 

IV. ANALOG QUANTUM SIMULATORS (AQS) 

In this section, we classify general properties of AQSs, 
discuss various proposals for such devices, and summa- 
rize the state-of-art knowledge concerning their control- 
lability, reliability, and efficiency. AQSs are experimental 
systems that are designed to mimic the quantum dynam- 
ics of interesting quantum many-body models, typically 
using "always on" interactions between particles that are 
augmented by fast local unitary control. While by defini- 
tion they operate in continuous time and thus the Trot- 
terization problems do not concern them, the standard 
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error-correction methods and fault tolerance cannot be 
applied. 

A. Universal Analog Quantum Simulators (UAQS) 



where, for simplicity, we shall term them AQSs. AQSs 
are experimental systems that can mimic continuous- 
time unitary Hamiltonian evolution for a family of mod- 
els of many-body physics. Their characteristics are as 
follows: 



Sometimes known as "Hamiltonian simulation," the 
goal of a UAQS is to transform a given Hamiltonian 
acting on a fixed Hilbert space into an arbitrary target 
Hamiltonian through a well-designed control sequence. 
While not conceived as a practical AQS device, the proto- 
col explores an abstract quantum-information-processing 
system capable of simulating unitary evolution for all (or 
at least all local) Hamiltonians. 

• Realizations: To our knowledge there are no 
concrete proposals for experimental realizations of 
UAQSs. 

• Controllability: While for UDQSs the issue is 
the access to the universal set of quantum gates, 
for UAQSs the question is what the necessary re- 
sources are (not necessarily quantum gates) that 
allow for the simulation of all Hamiltonian evo- 
lutions of interest. Universal control sets (as op- 
posed to universal digital logic gates) that gener- 
ate an arbitrary Hamiltonian evolution have been 
studied |136|, 1137] . Typically, such an approach us- 
ing "always on" interactions is associated with more 
limited control than is available in a universal dig- 
ital quantum computer. 

• Error correction: UAQSs do not allow for stan- 
dard error correction and fault tolerance. 

• Efficiency and Reliability: Diir et al. studied a 
hybrid construction of always-on interactions with 
stroboscopic digital control to achieve a universal 
Hamiltonian simulator via the Trotter construc- 
tion [138] . They found that decoherence and analog 
timing errors can make this inefficient for a Hamil- 
tonian simulator. Other issues concerning UAQSs 
are essentially the same as for non-universal AQSs, 
so we leave the discussion of them to the next sub- 
section. The only difference is that UAQSs, by def- 
inition, are capable of performing tests of robust- 
ness of the quantum simulations that we propose 
in the following section, i.e., tests based on adding 
disorder or noise in a controlled manner to the sim- 
ulated Hamiltonian. For non-universal AQSs such 
an addition requires additional resources. 

B. Non-Universal Analog Quantum Simulators 
(nUAQS), or simply AQS 

Non-universal AQSs constitute the most popular class 
of quantum simulators, but despite this fact, there is very 
little known about their reliability and efficiency. There- 
fore, we focus on them in the remainder of this paper. 



• Realizations: The most advanced experiments 
with AQSs have been with ultracold atoms in opti- 
cal lattices ISllS HT^fTOUTT] . The degree of quantum 
control is even better in ultracold-ion systems, but 
these are so far limited to few ions |21] - E51 [27] . 
The first step toward large-scale QSs with ions was, 
however, recently achieved |T39j. Recently there 
has also been substantial progress in investigations 
of other possible candidates for AQS, such as pho- 
tonic systems |31l [5^ I34n36i , photonic and polari- 
ton systems |33[ artificial lattices in solid 
state systems |44| . 

• Controllability: Most, if not all of the propos- 
als for and realizations of AQSs allow for at least 
partial controllability. The paradigm examples are 
AQSs employing ultracold atoms in optical lattices 
(for more details, see Chapter 4 of [3 ). Here, the 
typical controls involve optical lattice parameters 
(laser intensity, wavelength, etc.), lattice geometry, 
lattice dimensionality, temperature and other ther- 
modynamical control parameters, as well as atomic 
interaction strength and nature (van der Waals in- 
teractions are controlled via Feshbach resonances, 
while dipole interactions by the strength of the 
dipoles, lattice-site-potential shape, etc.). Further, 
tunneling can be laser assisted and can mimic arti- 
ficial Abelian or even non-Abelian gauge fields (cf. 

|T7]). Dipolar interactions may lead to non- 
standard terms in Hubbard models, such as occu- 
pation dependent tunneling |140| . or various effects 
involving higher orbitals (see, e.g., |141j ). 

• Error correction: AQSs do not allow for stan- 
dard error correction and fault tolerance. 

• Efficiency: The issues of reliability and efficiency 
are essential for the usefulness of any QS, and AQSs 
in particular. In the context of AQSs, however, 
there has been little analysis of these problems. 
Firm criteria on computational complexity and ef- 
ficiency for AQSs are in general difficult to address 
and have not yet been established. First of all, 
they require the knowledge of classical computa- 
tional complexity of the static or dynamical prop- 
erties of the considered quantum models. Unfortu- 
nately, in the realm of classical computation, there 
are few proofs that a given computational problem 
is outside the class P, or even if there is a clear de- 
lineation between certain complexity classes. In re- 
cent years, there has been considerable progress in 
understanding that the ground states of ID gapped 
systems can be efficiently simulated by classical 
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methods |142H144] . or that the quantum dynam- 
ics is in general computationally hard [Ml • If 
we can set the parameters of our AQS to a regime 
where efficient classical simulation is possible, we 
can assess the efficiency and reliability of the AQS 
in this case by direct comparison with classical sim- 
ulations (see below). However, there is no guaran- 
tee that such a calibration will hold in the truly 
interesting regimes of parameters, where efficient 
classical simulations are either impossible, or we 
do not know how to perform them. 

• Reliability: So far, there exists no perfect and 
rigorous way to assess the reliability of AQSs, but 
there are several complementary approaches. One 
approach is by cross validation of a variety of differ- 
ent physical systems (e.g., atoms in optical lattices, 
ions in traps, and superconductors) |74]. The hope 
is that since every platform has its own set of imper- 
fections, they will agree on the universal properties 
of the ideal quantum many-body model being sim- 
ulated. While it remains to be seen whether such 
universal features would emerge, this approach has 
a number of shortcomings. For example, there 
may be models that have only one known imple- 
mentation, or different implementations may suf- 
fer in the same way from imperfections, hence con- 
sistently exhibiting features associated with noise 
rather than with the ideal model. 

A more systematic approach is to validate results 
of a quantum simulator against analytical and nu- 
merical predictions in the regime of parameters 
where such comparison is possible. This was re- 
cently demonstrated in experiments with ultracold 
bosonic and fermionic atoms [6-8 ; amazingly, in 
one case numerical simulations helped to correct 
the expected experimental temperature by up to 
30%. Relying solely on validation from classical 
calculations, however, would restrict QSs to models 
in regimes where these efficient classical algorithms 
exist - that means contradicting our relevance and 
usefulness requirement, point (a) of our definition 
of a QS. In general, we desire to operate quantum 
simulators in regimes whose properties are difficult 
to deduce by classical methods, e.g., near or at the 
critical point of a QPT, or in genuine terra incog- 
nita regimes. In these regimes, however, many rel- 
evant models become hypersensitive to perturba- 
tions [I45J I146| , and even small levels of noise may 
spoil completely the results of the quantum simula- 
tion. Indeed, the capability of an analog quantum 
information processor whose dynamics is charac- 
terized by quantum chaos (i.e., well described by 
random matrix theory) can be severely impacted 
by imperfections |147II148] . More importantly, this 
also means that successfully validating a QS in a 
classically accessible regime does not give certainty 
about its robustness in regimes which are classically 



not accessible. 



C. Open-System Analog Quantum Simulators 
(OSAQS) 

Finally, let us mention Open-System Analog Quantum 
Simulators (OSAQSs). Similar to OSDQSs, OSAQSs are 
supposed to simulate dissipative dynamics described in 
the simplest situation by a continuous-time Markovian 
Lindblad master equation for the density matrix of a 
many-body system of interest. OSAQSs can be aimed at 
a simulation of interesting open-system dynamics, or at a 
designed dissipative dynamics toward a stationary state 
of interest. Many-body Lindblad master equations have 
been studied in the context of evaporative [ 149H155] , laser 
|156fil67] and sympathetic [168 -172J cooling of degener- 
ate atomic gases (see also |173j ). Recently, there has been 
a revival of interest in such systems in the context of pos- 
sibility of using them for preparation of interesting pure, 
highly-entangled states [126-128 . Open-system quan- 
tum simulators employing superconducting qubits may 
also give insight into exciton transport in photosynthetic 
complexes [174]. 

The experimentalist who realizes an OSAQS, in con- 
trast to an AQS has to have to his/her disposal some 
non-unitary, dissipative quantum mechanism: in a sense, 
all designed cooling or entropy-reduction methods are of 
this sort. 

• Realizations: All AQS systems can, in principle, 
be used as OSAQSs. 

• Controllability: OSAQSs are typically not uni- 
versal in a sense similar to OSDQSs; they allow 
neither for simulating arbitrary (local) Markovian 
dynamics, nor do they allow for preparation of ar- 
bitrary states. Nevertheless, in most of the propos- 
als |l26fll28] or experimental realizations they al- 
low for a wide control of parameters, which in turn 
allows for simulation of open-system (Markovian) 
evolutions for wide families of open systems. 

• Error correction: For OSAQSs, it is not guaran- 
teed that error correction and fault-tolerant com- 
puting is possible in the sense defined in the sub- 
section on OSDQSs. 

• Efficiency and reliability: All of the above 
discussion concerning AQSs applies also OSAQSs. 
But, again due to the purely dissipative nature of 
the process, this type of simulation has a certain 
amount of intrinsic robustness and built-in "error 
correction" - this is particularly clear for the OS- 
AQS of quantum kinetic Ising models or Kitaev's 
toric code |173j . Still, as in the case of OSDQSs, 
most of these general aspects concerning OSAQSs 
have not yet been investigated systematically. 
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V. ROBUSTNESS OF ANALOG QUANTUM 
SIMULATORS 

All of the above considerations clearly lead to the fun- 
damental question: Can we trust quantum simulators? 
From what we have said, the rigorous answer to this ques- 
tion is "no", yet in practice we do tend to trust them, at 
least to some extent. 

In order to gain more trust in the results of quan- 
tum simulations, it is thus extremely important to design 
novel tests and certification of reliability and validity of 
QSs. In this section, which constitutes some of the most 
important results of this paper, we propose such tests, 
which we call tests of robustness of quantum simulators. 
Our tests consist in checking the robustness of QSs with 
respect to addition of imperfections, such as static dis- 
order or dynamical noise. This would then allow to (i) 
judge how strong the reaction of the QS with respect to 
these perturbations is, and (ii) might even open possibil- 
ities to extrapolate interesting observables to the ideal, 
zero-disorder limit. Such tests can also be applied to 
DQSs, but are particularly suited for AQSs. For exam- 
ple, in an implementation with trapped ultracold atoms, 
disorder can be increased in a controlled manner [175] . 

In the following of this key issue article, we use the 
example of the quantum Ising chain to substantiate our 
discussion of the reliability of AQSs and the relationship 
to the complexity/efficiency of the simulation. We study 
how imperfections affect the results of a AQS simulating 
that model, where, for simplicity, we assume quenched 
disorder as the only possible imperfection. In the future, 
it will be in particular interesting to also investigate the 
effects of dynamical noise, and the decoherence and re- 
laxation that occurs due to coupling with an environment 
(see also Ref. [176|). 

The quantum Ising model, which we describe in 
detail in Section VI, is exactly solvable, which allows 
us to explore regions with universal behavior such as 
second-order QPTs. In Section VII, we show that the 
ground-state expectation values of certain local observ- 
ables appear fairly robust under disorder, while this 
need not be true for the global many-body state of the 
simulator. In particular, disorder can have a significant 
effect on relevant quantities that one could hope to 
extract from the simulator, such as critical points and 
exponents, or - if the system is described by a conformal 
field theory (CFT) - its central charge [177\. Finally, we 
briefiy address the relationship between robustness and 
complexity by studying the dynamics of different thermal 
states after a quench of the Hamiltonian (Section VIII). 
We show evidence that QSs appear to work better in 
regimes that are classically easier to solve or simulate 
(high-temperature states), thus hinting at a connection 
between the amount of quantum correlations and the 
robustness of a QS. 



VI. THE MODEL 

To illustrate the influence of disorder on an AQS, we 
study the transverse Ising model (TIM) 

H^-J2j,,a^a^-J2h^<7!, (1) 

where cr^'^ are the usual Pauli spin matrices and 
means sum over nearest neighbors. The system is sub- 
ject to quenched disorder in both the interaction and 
field terms. We denote the nearest-neighbor spin cou- 
pling and the transverse field by Jij = J(l -|- rSij) and 
hi = h(l -\- rrji), respectively, where dij and rji are inde- 
pendent random variables with a Gaussian distribution of 
mean zero and variance r. All details of our calculations 
are presented in the Appendix. 

The TIM, even under the presence of disorder, is effi- 
ciently solvable - by which we mean that the eigenstates 
and eigenenergies of the system can be found using a clas- 
sical computer, and that the cost of the algorithms (in 
time and hardware) is polynomial with the size (number 
of particles) of the system. The TIM, in particular, can 
be solved by using a Jordan- Wigner transformation to a 
system of non-interacting fermions - and the cost of solv- 
ing the non-interacting fermion system is the cost of di- 
agonalizing a matrix with rank equal to twice the number 
of spins in the chain [178J. This model is well studied (see 
for instance |179| ): for low fields, the ground state is a 
ferromagnet, while for large fields it is a paramagnet. At 
zero temperature and disorder, the system undergoes a 
QPT when the dimensionless control parameter X = h/J 
approaches the critical value, Ac = 1, i.e., when the field 
intensity equals the interaction strength. The infiuence 
of disorder can have dramatic effects on this phase dia- 
gram: imperfections can create new phases, or even de- 
stroy the ones we want to investigate. Indeed, in the 
TIM when the disorder strength is comparable with the 
interactions, the critical point disappears and is replaced 
by a so called Griffiths phase |180j . extending across a re- 
gion of size proportional to the disorder strength. Even 
more, in this Griffiths phase observables become non-self- 
averaging, i.e., fluctuations increase with system size, and 
hence dominate the thermodynamic limit. In this study, 
we consider small disorder strengths, which allows us to 
ignore the Griffiths phase, especially in finite-size sys- 
tems. Moreover, a state-of-the-art AQS can achieve very 
low levels of disorder, whence this is the experimentally 
relevant regime. Note that, while the TIM has been stud- 
ied extensively in the limit of very large disorder |.181J , 
there are few studies addressing directly the influence of 
small disorder on the universal properties near quantum 
phase transitions. However, it is well known that even 
small disorder can lead to novel quantum phases such as 
Bose fl82l or Fermi [183] glasses (see also the section on 
disorder in Ref. |184| and references therein). In the fol- 
lowing, we analyze the robustness of relevant observables 
to disorder in static and dynamic situations. 
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Figure 1. A: The correlation length ^ decreases with disor- 
der r, and its peak broadens (shown for a chain of 400 spins). 
The critical point (as extracted from a finite-size scaling of 
the energy gap A) moves to larger A with increasing disorder 
(black line). B: For a chain of 400 spins, we show the mean 
single- and two-site reduced simulator fidelities (/i and /2), 
and the total simulator fidelity (F) for a fixed disorder level 
r — 0.1. Local fidelities are more robust, which gives hope 
that local quantities can be reliable even if disorder deteri- 
orates the overall ground state. As expected, disorder has 
more severe effects close to the QPT. C: The central charge 
c (full circles, left axis), extracted from a fit to the part-chain 
entropy, and the critical exponent i/ (open circles, right axis), 
extracted from a collapse of the correlations in different chain 
lengths. Both change with disorder, which can lead to erro- 
neously assigning the QPT to an incorrect universality class. 
However, the change begins relatively smoothly at low levels 
of disorder. 

VII. RESULTS STATICS 

First, we investigate static properties of the AQS and 
their robustness to disorder (summarized in Fig.[lJ. We 
average all the analyzed static quantities over many re- 
alizations of disorder. 

One can evaluate the response of the AQS to disor- 
der using the simulator fidelity, which we define for pure 
states as the overlap between the state obtained with a 
perturbed simulator, |^'r(A)), and the ideal state |\I'o(A)), 

F(r,A) = |(vl/o(A)|vI/.(A))|. (2) 

Although we define the simulator fidelity for any pos- 
sible target state, we focus on the ground state. As 
Fig. [Tj3 shows, this overlap is considerably suppressed 
near the QPT, reaching values as low as 55% (for r = 0.1 



in a chain of L = 400 sites). When scaling to larger 
systems, F{r, A) will typically vanish exponentially fast, 
simply due to the exponential growth of the dimension of 
the Hilbert space (a kind of "orthogonality catastrophe"). 
In a universal quantum computation, the fidelity would 
have to be very close to 1 for the quantum computer to 
work fault-tolerantly. However, QSs have the advantage 
that we do not necessarily demand of the entire state to 
be robust. Often, it is enough if we can distinguish the 
relevant phases by measuring faithfully local observables 
(local in the quantum information sense that few sites 
are involved, although they may be physically far apart). 
Obviously, this is less demanding, yet very useful. 

To quantify the robustness of local observables, we in- 
vestigate the single-site and (nearest-neighbors) two-site 
simulator fidelity /i(r. A) and f2{r,X), respectively. As 
the one- and two-particle density matrices will generally 
be mixed when the overall pure many-body state is en- 
tangled, these are defined as the Uhlmann fidelity |185] 
between the single- or two-site reduced density matri- 
ces of the ideal state and the one at disorder strength r, 
/ = Tiy^yj^pry/pQ. It can be assumed that fidelities of 
the reduced system decrease more or less monotonically 
with the number of sites involved. As seen in Fig. [lj3, 
the reduced simulator fidelities are much more robust to 
disorder than the global one - near the phase transition, 
/2(r, A) decreases to approximately 0.998, and /i(r, A) re- 
mains above 0.999. This gives optimism that local quan- 
tities are robust enough to allow a faithful distinction 
between different quantum phases. 

One step beyond local properties of the ground state 
are the correlation lengths dictating the exponential de- 
cay of long-distance correlation functions. We investigate 
the correlation length extracted from the correlation 
function 

(3) 

where away from criticality C{i,j) oc exp(— |i — j | /^). 
Without disorder and for infinite systems, ^ diverges at 
the critical point, because criticality is the emergence of 
collective phenomena involving infinite degrees of free- 
dom at all length scales. In practice, we can only deal 
with finite systems so that we cannot observe real critical- 
ity but only smoothed out signatures of it, a phenomenon 
which one normally calls "pseudo-criticality." For exam- 
ple, the correlation length ^ is bounded by the system 
size. Still, its peak gives a reliable signature for the lo- 
cation of the critical point. As Fig. [iJV shows, however, 
disorder suppresses correlations and broadens the peak of 
^, thus making an extraction of the critical point much 
less reliable. 

Another criterion to locate the QPT is provided by 
the energy gap A between ground state and first ex- 
cited state. At criticality, the low-energy spectrum of 
the Hamiltonian is gapless in the thermodynamic limit. 
In finite systems, it presents non-vanishing gaps that de- 
crease in a systematic way with increasing system size. 
Due to this characteristic scaling of physical observables 
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as a function of system size in pseudo-critical systems, 
criticality can be detected by studying a sequence of fi- 
nite but increasingly large systems, a technique called 
finite-size scaling [186 . We describe this technique for 
the energy gap in the Appendix (see Fig. Al I, and show 
in Fig. [l|\ (black line) the location of the critical point 
extrapolated in this way. As can be seen, if one does not 
correct for disorder effects, one would locate the critical 
point at values of A that are too large. 

Perhaps of more fundamental interest than the ex- 
act location of a critical point is its universality class. 
All models within a given universality class give rise to 
the same collective behavior at large distances (typically 
large with respect to the lattice spacing), irrespective of 
their microscopic details |177| . Therefore, all relevant 
thermodynamic quantities for all models within a class 
are characterized by the same small set of critical expo- 
nents which describe the power-law decay of the corre- 
lation functions of local observables in the large-distance 
regime, a property that allows to differentiate among dif- 
ferent emerging collective behaviors. To investigate how 
robust the universal behavior is, we compute the critical 
exponent for the correlation length, i', from a collapse 
of the correlations (as explained in the Appendix, see 
Eq. ( |A10[ ) and Fig. [Alj ). As shown in Fig. [ip, already 
for a few percent of disorder, ly increases strongly from 
its ideal value 1. Therefore, if one simply neglects the 
influence of disorder, the extraction of critical exponents 
yields wrong results. 

If the QPT is described by a CFT (a specific subclass 
of one-dimensional critical systems), it is characterized 
by a central charge c. The central charge appears ubiq- 
uitously [187j . It, e.g., governs the temperature depen- 
dence of the free energy (Stefan-Boltzmann law) and the 
Casimir effect in finite geometries, but also the scaling 
of the entanglement entropy of sub-regions of the ground 
state of the corresponding quantum models. Models with 
different central charge have different emerging collective 
behavior. For example, models whose collective behavior 
is that of a free Majorana fermion (as in the disorder-free 
TIM) have central charge = 1/2, while models whose 
collective behavior is that of a free boson have central 
charge = 1. Strictly speaking, the TIM has an under- 
lying CFT only in the disorder- free case, but there have 
been efforts to extract an effective central charge also for 
the disordered model [181', for example, from the von 
Neumann entropy S of the reduced density matrix of a 
part of the chain of size I. At criticality, this entropy 
scales as c/6 log(i/7r sin(/7r/L)) + A [USHISS]. Figure[lp 
shows that disorder decreases the effective central charge. 
Hence, ignoring the effects of disorder would give com- 
pletely erroneous results, since even a small deviation of 
the central charge indicates a completely different univer- 
sality class. Note also that the decrease of c with disorder 
indicates the destruction of correlations by disorder. 

Fortunately, for all the extracted quantities (except the 
global simulator fidelity) , the levels of disorder for appre- 
ciable changes to occur are at least a few percent. If the 



AQS can be operated below such a value, its results seem 
to be robust, at least in this simple model system. State- 
of-the-art experiments are good enough to fulfill this re- 
quirement. In fact, in many experimental situations one 
hopes to reach levels of disorder or noise that are below 
a few percents, assuring the robustness of the AQS. Fre- 
quently, however, changing parameters from the regime 
where validation via classical simulation is possible to 
the regime of terra incognita might lead to uncontrolled 
disorder or noise. That is why checking sensitivity to dis- 
order and noise in those regimes where it can be checked 
is of great importance. 



VIII. RESULTS - DYNAMICS 

Efficient classical algorithms for computing static prop- 
erties of quantum systems are more developed than for 
computing dynamics (the difficulty arises mainly because 
entropy and correlations grow rapidly with simulated 
time). Therefore, one can assume that in the absence 
of disorder, a quantum simulation of dynamics can much 
more easily outperform classical computers. Indeed, in 
a recent experiment based on ultracold bosonic atoms, 
the controlled dynamics ran for longer times than present 
classical algorithms based on matrix product states could 
efficiently track . We thus turn to the issue of how 
disorder affects the reliability of quantum simulations of 
dynamics. As with statics, we investigate the behavior of 
the simulator fidelity, but now also as a function of time, 
initial state, and external driving. 

Typically, we expect that the simulator fidelity will 
decay with time, and eventually reach an asymptotic fi- 
nite value. The effect of disorder in both the decay rate 
and the asymptotic saturation value can, in general, be 
understood from established techniques such as Fermi 
Golden 's rule, and random matrices |191] . On the other 
hand, the effect of the initial state and the external driv- 
ing is known to be nontrivial and of particular interest 
for our purposes. For example, it is known that numeri- 
cal techniques such as the time-dependent density matrix 
renormalization group (tDMRG) can simulate efficiently 
the dynamics after a sudden quench of the field h, as 
long as the quench is restricted to a few sites on the 
chain. However, if the quench is global, it has been shown 
that the computational resources needed to keep a fixed 
amount of error grow exponentially with time |192l I193| . 
Generically, solving for the dynamics of a quantum many- 
body system is a hard problem for classical algorithms. 
Our model is special because it can be solved exactly 
for all cases, although it remains hard for the tDMRG 
algorithm. We use this to our advantage to study how 
this class of classical algorithms behaves when solving for 
quantum dynamics. 

We studied the behavior of the full simulator fidelity 
under the following driving. As initial state we prepare 
the ground state of the Hamiltonian for a given value of 
the external field. At time zero, the field is quenched 
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instantaneously to a larger strength, and the system is 
allowed to evolve. In panels B and C of Fig. [2] we com- 
pare the short- and long-time behavior of fidelity for the 
case of a global and a local (single-site) quench. The AQS 
keeps a high fidelity in the case of a local quench, while 
it performs poorly for the global quench, with fidelities 
reaching lows of 0.8 even for small systems of 50 spins. 
We also observe that the AQS performs worse when the 
quench crosses the critical point, as shown in Fig. |2p, 
where we fix the strength of the quench and vary the 
initial field value. 

The initial state can also have an effect on the effi- 
ciency of classical algorithms. Using the same setup with 
a global quench, but starting from a thermal initial state, 
tDMRG becomes efficient for high temperatures [192_ 
where the state and its correlations are almost classi- 
cal. However, it becomes exponentially inefficient with 
time for low-temperature initial states. For initial ther- 
mal states, we can still compute the dynamics exactly, 
although computationally it becomes too expensive to 
calculate the full many-body fidelity between the evolved 
states. In this case, therefore, we focus on the reduced 
simulator fidelity. For the regimes of disorder that we 
studied, we observe that the time dependent fidelity de- 
cays with a rate roughly proportional to the strength 
of the disorder squared (typical of a Fermi golden rule 
|191| 1. For this reason, we show in Fig. [2jV a rescaled 
form of the fidelity, (1 — /i)/r^, that exemplifies the typ- 
ical behavior for all disorder strengths, as a function of 
time and temperature of the initial state. 

As with the classical algorithms |192| . the AQS remains 
faithful when the state is almost classical (high temper- 
atures). The simulator fidelity decreases rapidly for low 
temperatures, although it saturates at a fairly high value. 
In terms of distinguishability, the values we find imply 
that a fair observer would have only a 4% chance of dis- 
tinguishing the 1-spin reduced state of the AQS from the 
ideal state. In the inset of the top panel we show the 
average asymptotic fidelity as a function of temperature 
of the initial state. Again, for low temperatures fidelity 
worsens, but saturates to a few percent. For high tem- 
peratures, it is simple to perform an expansion of the 
fidelity which shows that /i ~ 1 — . 



IX. DISCUSSION AND OUTLOOK 

A key issue for future investigation is the relationship 
between the robustness of an analog quantum simulator 
and its computational power. For the models we have 
considered here, the physically relevant correlation func- 
tions are robust for a reasonable degree of disorder. This 
suggests that such an AQS could perform well in a labo- 
ratory demonstration. But, the TIM that we considered 
here is simulatable on a classical computer. Is this con- 
nection between robustness and classical simulatability 
coincidental, or does it refiect a deeper relationship? 

Disorder reduces the correlation length of the spin 
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Figure 2. A: Evolution of the average reduced simulator 
fidelity as a function of the temperature of the initial state. 
The system is an Ising spin chain of length 50, the initial state 
is a thermal state at criticality (A = 1), and at time zero the 
field is suddenly quenched to A = 2. In the vertical axis 
we show the infidelity (one minus fidelity) normalized by the 
disorder strength r squared. For larger temperatures (where 
there are less correlations) the state is more robust. In the 
inset, we show the average asymptotic infidelity as a function 
of temperature. For large temperatures it decays as l/T"^. 
B: Evolution of the full simulator fidelity for an initial state 
equal to the ground state (zero temperature) at A = 0.75 after 
a sudden quench to A = 1.25. For a local quench in a single 
site, fidelity saturates rapidly at large values, but decreases 
strongly for a global quench. C: Asymptotic value of the total 
simulator fidelity as a function of the initial value of the field 
A, with a fixed quench strength of (5A = 0.25. The system 
is much less robust for global quenches and near criticality 
(A=l). 



chain. Because less-correlated quantum states can be 
described with fewer parameters, there is reason to sus- 
pect that certain aspects of weakly disordered quantum 
many-body systems could actually be easier to simulate 
on classical computers than their clean idealized versions. 
This happens, for example, in the realm of digital quan- 
tum computation, where a quantum circuit becomes clas- 
sically simulatable for noise levels above a certain level 
when quantum gates lose their entangling power |1941 - 
'196]. In the context of many-body physics, the success 
of DMRG in, e.g., ID spin chains, is rigorously related to 
the existence of efficient matrix-product-state represen- 
tations |197| . These take advantage of the small amount 
of quantum correlations in such systems, thus compress- 
ing the 0(exp(n)) parameters needed to describe a gen- 
eral n-particle state to 0{n) finite-dimensional matrices. 
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In higher-dimensional lattices, states which obey the so- 
called "Area law" [97,, where quantum correlations are 
smaller than in generic states, may still be amenable to a 
classical simulation using state-of-the-art techniques such 
as tensor networks |771 [HI [HI [Ml [HH], Density Func- 
tional Theories |86], or Quantum Monte Carlo |75| . 

We thus arrive at the fundamental question: Do the 
finite imperfections of an analog quantum simulator re- 
duce the correlations, and thus the number of parame- 
ters needed to describe the system, so as to render the 
device simulatable by classical means! We know that 
for noise above certain levels a digital quantum circuit 
is classically simulatable and for levels below a certain 
threshold it can be rendered fault tolerant. Is there an 
intermediate regime for which noise is too great to al- 
low fault-tolerant universal quantum computation, but 
small enough that an AQS accesses physics beyond clas- 
sical simulation? The existence of an intermediate regime 
would imply that there exists a whole class of problems 
outside P that we can access in the near future, even 
without a fully functioning quantum computer. 

The results we present here, in particular those for dy- 
namics, are an initial attempt - albeit in a trivial model - 
at understanding the above problem. We can see how an 
analog quantum simulator works well when a classical so- 
lution is efficient, and worsens (but only in a limited way) 
when the problem becomes classically hard to simulate. 
Even though the underlying model is actually solvable, 
this may be positive evidence for the existence of an in- 
termediate regime of noise, and the efficiency of AQSs in 
more complex situations. 

Our main discussion focused on AQSs, but similar is- 
sues pertain to DQSs. Since to date there exists no 
known way to fault-tolerantly error-correct AQSs, there 
is a natural tendency to explore the advantages of DQSs, 
where error correction is possible. The above discussion 
shows, however, that a digital implementation of a quan- 
tum simulation does not, in itself, guarantee an efficient 
and more powerful simulation than one that is carried 
out classically. As in any quantum algorithm, initial- 
ization, evolution of the state, and measurement must 
be performed efficiently, i.e., with a polynomial use of 
physical resources (space and time). Digital quantum 
simulation is no exception. Indeed, as discussed above, a 
fault-tolerant implementation of the standard approach 
based on the Trotter expansion [M^ comes at the cost of 
an overhead in the number of gates and time required 
that grows exponentially with the degree of precision re- 
quired [ 641 1120] . If we can guarantee the reliability of 
analog quantum simulators while avoiding such exponen- 
tial costs, many open problems from all areas of physics 
could suddenly come into the reach of being solved. 

Finally, we can turn the problem of quantum simu- 
lation on its head and ask, what does Nature do? For 
any real material, like a high-Tc cuprate, has imperfec- 
tions. Does Nature access highly correlated states that 
cannot be efficiently simulated on a classical computer? 
Certainly, in some cases we believe it does, as for exam- 



ple in high-Tc superconductors [SS] or in certain ground 
states of frustrated quantum antiferromagnets which are 
believed to carry topological order |SS]. If noise is low 
enough, does Nature protect quantum correlations to a 
degree that classical methods cannot efficiently represent 
the physically interesting quantities? And, can we exploit 
this capability with a quantum simulator? If Nature does 
it, we should take advantage of it! 
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Appendix 

Quadratic fermionic systems — The transverse field 
Ising model, Eq. (1) of the main text, even with disor- 
der, can be solved by casting it into the form of non- 
interacting fermionic particles using the Jordan- Wigner 
transformation, 



m— 1 



m— 1 

2c]c, 



(Ala) 

(Alb) 
(Ale) 



The Cj, c' obey on commutation relations. This trans- 



formation leads to 



H = Y. [^A^J'^J + {4B^3C\ -f h.C. 



2 X! ' 

3 



where 



Hamiltonian ( A2 1 can be diagonalized to 



H 



N 

^ J^kvivk + Eq , 
fe=l 



(A2) 

(A3a) 
(A3b) 



(A4) 



where A — ^{A — B)'^'' is diagonal. A, <I>, and 
can be obtained from the singular-value decompo- 
sition of Z = A — B. The normal modes are rjk = 

J2f=i idk^jCj + hk,jcfj , where g = ($ + *) /2, and h = 
((f> — 5*) /2. From this, we can compute the relevant 
ground-state properties. 
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Ground- state fidelity and correlations — From the 
normal modes obtained in the diagonalization of the pre- 
vious section, we can compute the observables we are 
interested in: the simulator fidelity F (the overlap to the 
disorder-free ground state), reduced simulator fidelities, 
the energy gap, and the ZZ-correlations. 

In general, the overlap between the ground states of 
two realizations Z and Z is [200j 



F{Z,Z] = Wdet 



1 + T-^T 



(A5) 



with T = ($~^A$)~ Z. We define the simulator fidehty 
F as the overlap at fixed A between the ideal, disorder- 
free state and the state at disorder strength r, 



F{r,X) = F{Z{X)r,Z{X)o) 



(A6) 



This is a global quantity, but one can expect that local 
observables are less affected by disorder. A measure for 
the change of local quantitites is the single-site simulator 
fidelity 



h {r, A) = ti-^\f^)pi'\x)\f^) . 



(A7) 



(i) . 

where pr = ti j^ip is the reduced density matrix of 

site i under disorder r, and Pq ^ is the equivalent in the 
disorder-free case. The single-site reduced density matrix 
is completely determined by the expectation values of trf , 
p — x^y^z, since one can expand p^*) = | X);^ ('^f ) 
Here, the sum runs over cr'^, p = x,y,z, and cr*^'^-' — 1. 
We also analyse the two-site simulator fidelity 

/2 = E tr/y^p^p(-+^) (A) vS^-^+^) (A) , (A8) 



for nearest neighbors. Here, pi^'^^^^ = tr^^^i ^+ip is the 
reduced density matrix of sites (i, « + 1) under disorder r, 

and Po is the equivalent in the disorder-free case. We 
compute all considered static quantities as the mean over 
a large number of disorder realizations; for the fidelities 
F,fi, and /2 displayed in Fig. IB of the main text, we 
used 5000 realizations at chain length L ~ 400. 

The correlations, finally, can be computed using the 
fact that the ground state j'l') of Eq. (A4| is the vacuum 
of the normal modes (i.e., rjk j^f) = 0, Vfc). For example, 
for the ZZ-correlations this yields 



Ci^,J) ^ i^M'^^i'^l^r) - {^M'^\^r){^M'^\^r) 

= 4(*r|clc,c]cj|*,.) - 4(*^|4c,|*^)(5'^|c]cj|^'0 
= 4 (h^h)^^ (g^g)^^ - 4 (h^g)^^ {g^ . (A9) 

Away from criticality, the correlations decay as C{i,j) oc 
exp(— \i — j \ /£,) with correlation length ^. In Fig. lA of 
the main text, we display ^ extracted from fits to part 



of the wings of C{i,j) (for L = 400 and 10000 disorder 
realizations) . 

Without disorder, 



C{^,J)L'^' ^fi\t-j\/L) 



(AlO) 



for some universal function / |177| . Hence, one can ex- 
tract the critical exponent for the correlation length v 
from a data collapse of the correlations. In Fig. IC of 
the main text, we show the erroneous values for i', ex- 



tracted from Eq. (AlO I if one naively neglects that this 



relationship is no longer true in the presence of disor- 
der. Figure |A1| shows the best collapse achieved with 
Eq. (AlO) for disorder levels r = and 0.2. The value 
of 1/ for the best collapse increases with disorder. Hence, 
using Eq. (AlO I on a disordered AQS yields a too large 
critical exponent, compared to the ideal model. More- 
over, the qualitity of the collapse worsens with increas- 
ing disorder, demonstrating that a naive application of 
Eq. (AlO I is unjustified if disorder is large. For this anal- 
ysis, we used L = 100 to 190 in steps of 10 with 10^ 
disorder realizations, L = 200, 250, and 300 with 5 x 10^ 
realizations, and L = 350 and 400 with lO'^ realizations. 

The correlations are intrinsically connected to the en- 
ergy gap A(L), since a gapped system necessarily has ex- 
ponentially decaying correlations [201] . Via a finite-size 
scaling, the gaps at finite systems also allows to extract 
the location of the QPT in the infinite system, as seen in 
Fig. Al There, we plot curves 1/(L^A(L)) for closeby 
chain lengths L, where C is the dynamical critical expo- 
nent, which for the disorder-free case equals 1. These 
cross at a series of pseudo-critical points which with in- 
creasing L tends rapidly to the critical point of the ther- 
modynamic limit [202] . Assuming that it does not change 
much for small disorder, we identify as an approximation 
the critical point with the mean of the crossing points 
curves 1/(LA(L)) for L = 100,150,200,300,350,400 
with 10000 realizations of disorder each. As displayed 
in Fig. |A1| the crossing point moves to larger values of 
A with increasing disorder (see also Fig. lA of the main 
text). This means that applying this analysis to a real- 
world AQS without correcting for disorder can yield erro- 
neous results (compared to the ideal model) for the loca- 
tion of the QPT. Moreover, the crossing point becomes 
less well defined with increasing disorder showing that 
this analysis should be corrected for the presence of dis- 
order. Finally, at large disorder, a second crossing point 
appears at lower A. The two crossing points open up to a 
V-like structure with increasing disorder. This could be 
interpreted as an indication of the Griffiths phase (the 
crossing points are qualitatively consistent with the ex- 
tent of the Griffiths phase found in [203 ). For a more 
quantitative analysis, however, one would need to ac- 
count for a change of C with increasing disorder. 

If a given model can be described by a conformal field 
theory, its universal critical behavior (including all criti- 
cal exponents) is completely defined by a single number, 
the central charge c. To extract it, we compute the von 
Neumann entropy S of the reduced density matrix of a 
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Figure Al. Upper row: For low disorder (r = 0, left), 
the curves C{i,j)LF''^ for different L plotted as functions of 
\i — j\ / L collapse perfectly for the physically correct value of 
the critical exponent v = 1 (dots of different color correspond 
to different L). For increasing disorder {r — 0.2, right), the 
collapse worsens, and the best collapse is obtained for some 
v > 1. Lower row: For low disorder (r — 0, left), the 
curves 1/(LA(L)) cross perfectly at the location of the critical 
point, A = 1. With increasing disorder (r = 0.2, right), the 
crossing point moves to larger values of A and becomes less 
well defined. Also, at large disorder there appears a second 
crossing point below A = 1. The insets show zooms on the 
crossing points. The findings of both rows of figures mean 
that, ignoring the effects of disorder in an imperfect AQS 
yields, compared to the ideal model, too large values for the 
critical exponent as well as the location of the critical point. 

part of the chain of size Z for L = 300 with 10000 disorder 
realizations. For systems with open boundary conditions, 
a fit to c/Q\og{L/TT am{lTT/L))+A (excluding small values 
of /) yields the effective central charge c |188fll90] . The 
disorder-free value is c = 0.5. Increasing disorder sup- 
presses this, indicating the decrease of entanglement in 
the system (see Fig. IC of the main text). Again, apply- 
ing the analysis that is correct in the disorder-free case 
(where the system is indeed described by a conformal 
field theory) without adjustments to the disordered sys- 
tem, yields results which deviate from the desired ideal 
case. 

Time dependent fidelities — For time evolution, we 
distinguish between the zero and finite temperature fi- 
delities, although the underlying technique is the same. 
We start by rewriting the fermionic Hamiltonian above 



as 

if = ivj/t .H.v]/, (All) 

where = (c| , cjy, Ci, cjv) is a 2A'' length vec- 
tor composed of all creation and annihilation operators 
present in H, and H = ^ (g) CTz -|- (g) cr^ is a 2N x 2N 
matrix with complex coefficients. 

For computing fidelities, we use the convenient Lev- 
itov's formula |204) I205| . which relates traces of oper- 
ators in the Hilbert space of the fermions to determi- 
nants of much smaller matrices (like H). For example, 
let P = vj/t . p . v[/ and Q = • Q • * be two operators 
in the space of fermions, with P and Q complex valued 
2N X 2N matrices. Then, 

Tr (e^e'^) = det (l -I- e^e^) . (A12) 

Similar formulas hold for more or less operators. 

In the zero-temperature case, when the initial state 
remains pure after evolution, the fidelity takes the form 
of an overlap 

F=\{Mt)\Mt))l (A13) 

where \tpo{t)) = e^'^o* j^^^^ jg ^^j^g initial state evolved 
with the target Hamiltonian of the simulation, Hq, and 
IV'r(O) = g-^^r-t jg ^]-^g game state evolved with an 
imperfect Hamiltonian H,- = Hq + rV . Rewriting the 
fidelity, 

F ^Ti poe'"^^e''"°\ (A14) 

with po — IV'o) (V'oIj we can use Levitov's formula and 
obtain 

F = dot (1 - Go + Goe'^-'e^*""*) , (A15) 

with Go = ("001 G|V'o)j and G the correlation matrix of 
the original fermionic operators, Gi j = 

If the initial state is not pure, but a thermal state, 
the state remains mixed even if the evolution is unitary. 
In this case, we cannot compute the fidelity for the full 
many-body state, but only the fidelity of the reduced den- 
sity matrix for a few spins. For this we must evaluate the 
correlation functions of the Pauli operators at different 
sites of the chain. For the case of a single spin, the sym- 
metry of the system ensures that at all times the reduced 
density matrix can be written as p = (1+ < a\. > a\.)/2. 
Since < cr* >= c|^Ci, we only need to compute the evolu- 
tion of the diagonal terms in the G correlation matrix. 
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